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ABSTRACT 

Infrared carbon stars without visible counterparts are generally known as ex- 
treme carbon stars. We have selected a subset of these stars with absorption 
features in the 10-13 /im range, which has been tentatively attributed to silicon 
carbide (SiC). We add three new objects meeting these criterion to the seven 
previously known, bringing our total sample to ten sources. We also present 
the result of radiative transfer modeling for these stars, comparing these results 
to those of previous studies. In order to constrain model parameters, we use 
published mass-loss rates, expansion velocities and theoretical dust condensa- 
tion models to determine the dust condensation temperature. These show that 
the inner dust temperatures of the dust shells for these sources are significantly 
higher than previously assumed. This also implies that the dominant dust species 
should be graphite instead of amorphous carbon. In combination with the higher 
condensation temperature we show that this results in a much higher acceleration 
of the dust grains than would be expected from previous work. Our model results 
suggest that the very optically thick stage of evolution does not coincide with the 
timescales for the superwind, but rather, that this is a very short-lived phase. 
Additionally, we compare model and observational parameters in an attempt to 
find any correlations. Finally, we show that the spectrum of one source, IRAS 
17534—3030, strongly implies that the 10-13 /xm feature is due to a solid state 
rather than a molecular species. 



Subject headings: infrared: stars — stars: carbon — circumstellar matter — dust 
— stars: AGB and post-AGB 
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1. Introduction 

1.1. Physical evolution of intermediate mass stars and their circumstellar 

shells 



St ars between about 1 an d 8 Mq will eventually evolve up the Asymptotic Giant Branch 
(AGB; llben &: Renzinilll983l ). Because of instabilities in their i nterior, AGB stars pulsat e 



and throw off large amounts of mass from their surface (e.g., iVassiliadis &: WoodI Il993l ). 
This intensive mass loss produces a circumstellar shell of dust and neutral gas. Once the 
AGB star has exhausted its outer envelope, the AGB phase ends. At this point, the mass 
loss virtually stops and the circumstellar shell begins to drift away from the star. At the 
same time, the central star begins to shrink and heat up from ~3000 K until it is hot 
enough to ionize the surrounding gas, at which point the object becomes a planetary nebula 
(PN). The short-lived post-AGB phase, as the star evolves toward to the PN phase, is also 
known as the proto- or pre-planetary nebula (PPN) phase. During the ascent of the AGB 



the velocity of the outflowin g mass appears to be fairly constant (e.g.. iHuggins et al.lll988 
Fong. Meixner fc Shahll2003l ). Therefore the dust furthest from the star represents the oldest 
mass loss, while material closer to the star represents more recent mass loss. Towards the end 
of the AGB phase t he increasing impact of th e thermal pulse cycles leads to an increasing 
mass-loss rate (e.g., IVassiliadis fc WoodI Il993l : IVillaver et al.l l2002al ). Such an increase in 



mass- loss rate (dub bed the superwind) is necessary to explain the densities seen in typical 



PNe ( lRenzinilll98ll ). Since the invocation of the superwi nd, many observations of AGB stars 



and p ost-AGB stars have supported this hypothesis (e.g.. lKnapp fc Morrislll985l : IWood et al. 
1992[ ). 



1.2. Chemical evolution of AGB stars and their circumstellar shells 

The chemical composition of the atmospheres of AGB stars is expected to change as 
these stars evolve, due to convective dredge up of carbon produced in the He-burning shell. 
The amount of carbon relative to oxygen (the C/0 ratio) is critical in determining which 
types of dust and molecules are present around an AGB star. The formation of extremely 
stable CO molecules will consume whichever of the two elements is less abundant, leaving 
only the more abundant element available for dust formation. Stars start their lives with 
the cosmic C/0 ratio of ?a0.4 and are therefore oxygen-rich. In about a third of AGB stars, 
enough carbon will be dredged up to make C/0 > 1 and therefore carbon will dominate 
the chemistry around these stars, known as carbon stars. Carbon stars are expected to 
have circumstellar shells dominated by amorphous or graphitic carbon grains, which do not 
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have diagnostic infrared features. Another component of the dust shell around carbon stars, 
silicon carbide (SiC), does have an infrared spectral feature at ^11/im and therefore has 
been of great interest to researchers s eeking to understand the evolution of the dust shells 



and infrared features of carbon stars (Baron et al. IQStI: Chan fc 



I995I: ISoeck et aljll997l: Isioan et al.l[l998l: ISpeck et al 



2005 



200 



<Cwok 



19901: 



Goebel et al. 



Thompson et al 



2006h. 



1.3. Extreme Carbon Stars 



As carbon stars evolve, mass loss is expected to increase. Consequently, their circum- 
stellar sh ells become progressivel y more optically thick, and eventually the central star is 
obscured. IVolk et al.l (119921 . 12OOOI ) christened su ch stars "extreme ca rbon stars" . These stars 
have al so been dubbed "in frared carbon stars" (IGroenewegeru Il994j ) . and "very cold carbon 
stars" (lOmont et al.l 119931 ). Extreme carbon stars are expected represent that small subset 
of carbon-rich AGB stars which are in the superwind phase, just prior to leaving the AGE. 
Because the superwind phase is short-lived compared to the AGB phase the number of ex- 
treme carbon stars is intrinsically small. Consequently, few of t hese objects are known. At 
present there are ~30 known extreme ca rbon stars in the Gala xy (jVolk et al.lll992l ) compared 
to ~30,000 known visible carbon stars (ISkrutskie et al.ll200ll ). 



van der Veen fc HabingI (119881 ) attempted to define a way to distinguish between oxygen- 



rich and carbon-rich AGB stars using IRAS color-color space, which was divided into sub- 
sections according to the pr operties of the dusty shells a re these star s (see Table 1 of 
van der Veen fc Habing] Il988l ) . This was further refined by lOmont et al.l (119931 ) who iden- 
tified a population of very cold carbon stars using HCN and CO observations, and showed 
that the regions originally designated as extremely dusty 0-rich AGB stars also contain a 
significant fraction of C-rich stars. 



The refinement of the lvan der Veen fc HabingI (119881 ) color-color diagram b ylOmont et a 



(I19931) defined subdivisions of the seven zones in color-color space (see Fig. 1 in lOmont et al. 



I993I ). Cool carbon stars with high mass-loss rates (and little or no SiC emission) fall into 
regions III and IV, which had previously been assumed to define OH-IR stars (i.e. the 
oxygen-rich counterparts to extreme carbon stars). The numbered regions have been sub- 
divided into smaller regions denoted by Illal, IIIa2, Illbl etc. A subset of the color-color 
space, covering parts of regions Illalc, IIIb2, IIIb2 and VIb is reproduced in Fig. [T] and 
includes our sample stars. 
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1.4. SiC absorption features 



SiC has long been pre dicted to be present in carbon star circumstellar sliells, beginning 
with condensation theories (jFriedmanlll969l : lGilmarull969[) and contin uing with the prediction 
of a characteristic SiC ~ll/im spectral feature (iGilra fc Codelll97ll ) and th en the observa- 
tional discovery of an ^ llp,m emission feature in many carbon star spectra fjHackwelllll972 



Treffers fc Cohenlll974j ). The effect of the evolving dust shell density structure on observed 
fe atures, and particu larly on the ~ll/im feature, have been discussed extensively (see review 



m 



Speck et al.ll2005l . and references therein). As the optical depth of the dust shell increases, 
self-absorption will diminish the ~ll/im feature and it will eventually be seen in net absorp- 
tion. SiC self-absorption was found t o be important in producing accurate radiative transfer 



models of extreme carbon stars (e.g. IVolk et al.lll992l ). even though this previous work did 



not recognize SiC absorption features. These absorption features are rare and have mostly 
been ignored in discussions of evolutionary sequences in carbon star spectra. In fact the 
rarity of such absorption fea tures led to the hypothesis that SiC bec omes coated in carbon 
at high optical depths (e.g. Baron et al. 1987 : Chan &: Kwokl E99ol ) . However, meteoritic 
data and theoretical models do not support this hypothesis (see § 11.61 and § 13.4.11) . 



A few extreme carbon stars have been shown to have an absorption feature at ~ 11/xm 
which has been tentatively attributed to SiC. This feature was discovered in the "prototype" 



extre me carbon star AFGL 306 8 (hereafter referre d to as IRAS 23166+1655; IJones et al. 



19781 ). and was re-examined bv ISpeck et al.l (119971). which a lso identified three additional 
extreme carbon stars with this feature. Clement et al. J2003h examined the absorption fea- 
tures of two of these extreme carbon stars (IRAS 23166+1655 and IRAS 02408+5458), and 
showed that their ll/xm absorption features are consistent with /S-SicJll nanoparticles. The 
broad absorption features of IRAS 19548+3035 and IRAS 21318+5631 (also discovered by 
Speck et al.l 119971 ) we re attributed to SiC abso rption with an interstellar silicate absorption 
contribution (see also Groenewegen et allll996l ). This will be discussed further in § 14.61 The 
ab sorption features in th e spectra of IRAS 19548+3035 and IRAS 21318+5631 were revisited 
by lClement et al.l (120051 ) who sug gested Si,sN4 grains as the carrier. H owever, this hypothesis 
has been shown to be erroneous (jPitman. Speck &: Hofmeisterll2006l ). 



The failure of the Si3N4 hypothesis led ISpeck et al.l (120051 ) to suggest that amorphous 
SiC grains may be able to account for the breadth, structure and barycentric position of 
the observed broad 10-13/im feature in IRAS 19548+3035 and IRAS 21318+5631. However, 
the dearth of amorphous presolar SiC grains seems to preclude this hypothesis (see § II. 6p . 



^Si li con carbide ex i sts in many (>70) di f ferent crystal structures, known as polytypes. See ISpeck et al 
|l997l ): lOaulton etaP |2003h : IPitman et aD |2008h for a discussion of the polytypes of SiC. 
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An alternative explanation for this feature is molecular line absorption, however, currently 



available line lists are not sufficient to properly assess this hypothesis (see lSpeck et al.ll2006 



and references therein). One molecular candidate which has trans itions in the correct wave- 



length range is C3 (e.g. IZijlstra et al.ll2006l : |j0rgensen et al.ll2000l ). but the line lists are not 



readily available. Furthermore, C3 is expected to be photospheric, rather than circumstellar, 
which probably precl udes its detection in op tically obscured stars. Moreover, the theoretical 
spectrum of C3 from I J0rgensen et al.l (I2OOOI ) shows a strong absorption close to the ~ 5/im 
CO line, which is stronger than the ~ 11/xm feature. As will be seen in § 14.41 the spectrum 
of IRAS 17534— 3030 does not show the 5/im absorption band and provides evidence that 
the observed absorption feature is not molecular in origin. 

Though previous research has included the effects of SiC self-absorption (shown to be 
crucial to produce accurate models Volk et al. 1992 : Speck et al. 1997 . 20051 ). no work has 
been done to directly fit the apparent SiC absorption feature in radiative transfer models of 
extreme carbon stars. 



1.5. Previous Radiative Transfer Models of Extreme Carbon Stars 



Volk et al.l (119921 ) perforn ied radiative transfer mo deling in order to match the Infrared 
Astronomical Satellite {IRAS; iNeugebauer et al.l 1 19841 ) Low Resolution Spectrometer (LRS) 
data for several of extreme carbon stars. They determined that the exact star temperature 
entered into the model was not importa nt for the emerging spectra due to the very thick dust 
shells around extreme carbon stars (c.f. iDePew et al.ll2006l : ISpeck et al.ll2000l ). Their models 
used a fixed composition ( a mixt ure of graphite and SiC), and a fixed dust condensation 
temperature. iGroenewegenI (119941 ) also performed radiative transfer modeling on a larger set 
of extreme carbon stars, but these models varied the dust c ondensation te r npera ture. Again 
this was based on IRAS LRS data. Following up on this, iGroenewegen (|l995l) modeled a 
large sample of carbon stars using amorphous carbon optical constants (iRouleau &: Martin 
19911 ). and assumed low dust condensation temperature in a fairly narrow range (650-900 K 
for the extreme carbon stars). Consequently the inner dust radius is larger than expected. 
Moreover the resulting models all have relatively low optical depths (rn.s^m < 2). The 
optical depth for their IRAS 23166+1655 mod el was found to be T ^ ^,:^f,m < 1, even though 



this star has an absorption feature at ll/im. IGroenewegen et al.l (119981 ) remodeled these 
stars, again assuming relatively low dust condensation temperatures, with similar results. 



Finally, IVolk et al.l (I2OOOI) used th e improved spectral resolution of the the Infrared Space 
Obse rvatory (/^O messler et al.lll996l ) Short Wavelength Spectrometer (SWS:lde Graauw et al, 
19961 ) to examine five extreme carbon stars. In their modeling study, IVolk et al.l (I2OOOI ) al- 
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lowed the optical depth and radial dust density distribution to vary; the resulting optical 
depths were relatively high (1.4-4.5 at 11.3/im), and the density of the dust shell was found 
to increase rapidly towards the center. This increase was interpreted as evidence of an in- 
creasing mass-loss rate within the last few thousand years, consistent with the identification 
of extreme carbon stars as the final stage of AGB star evolution. While these models did 
include SiC opacity data, the 11/xm absorption feature was not recognized and consequently 
no attempt was made to fit this feature in these models. 

A summary of the parameters of previous models for extreme carbon stars with ll/xm 
absorption features in our sample can be found in Table [H Interestingly, all previous models 
assume relatively low inner dust temperatures. This will be discussed further in § 13.4.11 
Furthermore, the modeled dust density distributions suggest a relatively slow increase in 
mass-loss rates (l/r^^, where x ~ 2.25 — 3.0). 

In reality, dust shells are expected to have heterogeneities and anisotropics in their den- 
sity structure as a res ult of pulsati on-driven dust formation and the ensuing hydrodynamic 
turbulent effects (e.g. IWoitkell2006l ). These dust formation models suggest that carbon star 
mass-loss is expect ed to be modu lated on several timescales, especially that of the pulsation 
cycle. Furthermore IWoitkd (120061 ) has suggested that the dynamics in the dust-forming zones 
around carbon stars lead to inhomogeneous dust formation, producing fine scale structure 
in the density of/, the dust envelope. In ad dition, while puls ation shocks are predicted to 



have a strong effect on local conditions (e.g. ICherchnefj 120061) . this is not reflected in tem- 



poral changes in the IR spectra of carbon stars (jCorman et al.ll2008l ). As will be seen in 
§ m the spatial scale of the heterogeneities is small and the timescale for pulsations is short 
compared to the timescales associated with even the thinnest dust shells. Moreover, the 
i nhomogeneities are e xpected to be wiped out over time by the hydrodynamic interactions 
( IVillaver et al.ll2002a| ]bl). Consequently, we do not consider these small scale structures in 
our models. 



1.6. Meteoritic Evidence 



The isotopic compositions of certain grains found in primitive meteorites indicate that 
they originated outside the solar system and are thus dubbed "presolar" . Dust grains from 
AGB stars are found virtually unaltered in these meteorite s, demonstrating that the se grains 
become part of the next generation of stars and planets (IClayton &: Nittlerl |2004J . and ref- 
erences therein). The precise physical characteristics of these meteoritic dust grains (e.g. 
sizes, crystal structures, compositions) can be used to help constrain the nature of the dust 
we see in our astronomical observations. 
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1.6.1. Presolar silicon carbide 



Silicon carbide was the first pres olar grain to be found in meteorites (IBernatowicz et al. 



19871 ) and remains the best studied (IBernatowicz et al.ll2006l . and references therein). The 
most important findings of this work are (1) that most (~ 99%) of the SiC presolar grains 
were formed around carbon stars; (2) of the AGE SiC grains, ^95% appear to originate 
around low-mass carbon stars (<3Mq), based on nucleosynthesis models of isotopic com- 
positions; (3) that all the SiC grains are crystalline (not amorphous); (4) that nearly all 
80%) are of the cubic /3 polytype, with the remainder comprising the lower temperature 
2H polytype; (5) that with one exception, SiC grains have not been found in the cores of 
carbon presolar grains (unlike other carbides: TiC, ZrC, and MoC); and (6) that the grain 
size distribution includes both very small and very large grains (1.5 nm — > 26 /xm), with most 
grains in the O.l-l/^m range. Single-crystal grains can exceed 20/im in size. Observations 
of the 11/im feature have been compared with laboratory spectra of various forms of SiC, 
and after some false starts it has now been attributed to /j-SiC, matchin g the information 



retrieved from meteoritic samples (jSpeck et al.lll999l : IClement et al.ll2003l ). However, there 



are still some discrepancies between observational and meteoritic evidence (most notably 
related to grain size). 



Prombo et al.l (119931 ) found a correlation between grain size and the concentration of 



s-process elements in SiC grains taken from the M urchison meteorites. The Indarch me- 



teorite presolar SiC grains yielded similar results (j Jennings et al.l |2002| ). In both cases 



the smaller grains have higher relative abundances of s-process elements. This observation 



(Laeadec et al. 


2007 




i008 


). 


dredge-up JSpeck et al 


2005) 



20081) . Alternatively, it may refiect an evolution in grain-size with 



1.6.2. Presolar "graphite" 
In addit ion to SiC presolar grain s, carbon grains are also relatively abundant and well 



studied (see IBernatowicz et al.l l2006l . and references therein). Presolar carbon grains are 
usually referred to as "graphite" grains, but their structures are more complex than this 
name infers. Presolar graphite is found in two types of spherules classified according to 
their external morphologies as "onion-like" and "caulifiower-like" . In general the graphite 
spherules follow a similar size distribution to the SiC grains. However, the high-density 
grains (p ~ 2.15 — 2.20 gcm^) associated with AGB stars have a mean size of 2/im. In 
addition, the AGB presolar graphite spherules span a larger range of isotopic compositions 
than the SiC grains, possibly suggesting that they form at a wider range of times during the 
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AGB phase. 

While the presolar SiC grains tend to be single crystals, the graphite grains regularly 
contain carbide grains. These carbides are enriched in s-process elements, indicative of 
formation around late-stage AGB stars. Many of the "onion-like" graphite grains have a core 
mantle structure in which the core contains disordered agglomerations of graphen^ sheets 
and PAHi-like products, while the mantle is composed of well-ordered graphitic concentric 
shells. The graphene particles have a typical size of 3-4nm. The "cauliflower-like" graphite 
grains also have a concentric shell structure, but it is less well ordered, and is composed 
primarily of the disordered graphene. Whether "onion" and "cauliflower" graphites are 
formed in the same outflows is not known. Both types of grain contain the refractory 
carbides and both span the same range of isotopic compositions. Whether the "onion" or 
"cauliflower" grains are more representative of grains in the outflows of extreme carbon 
stars is not known. However, even the most disordered "cauliflowers" or "onion" -cores are 
still closer to graphite than glassy carbon in structure. The least ordered grains are still 
considered to be agglommerations of nano-crystalline grains, rather than truly amorphous 
(pers. comm. K. Croat). 



1.6.3. Other presolar carbides 



As discussed in § 11. 6. H and § 11.6.21 refractory carbides are found inside "graphite" grains 
but not in SiC grains. Furthermore, SiC is not one of the carbides found in "graphite" grains. 
The refractory carbides (TiC, ZrC, MoC and RuC) provide more constraints on the dust 
formation processes around carbon stars. In particular, the formation of "graphite" spherules 
with TiC nuclei limits the range of C/0 ratios in which these grains could form to 1 ^C/0 
^ 1.2. Meanwhile, the ZrC can form nuclei at higher C/0, but the value still needs to be less 
than two. This is consistent with the measured C/0 ratios of Galac tic carbon stars, which 
have an average of 1.15 and a maximum of 1.8 (ILambert et al.lll986l ) 



^Graphene is basically a single sheet of graphitic material. If it is disordered, there are some heptagons 
and pentagons in place of the regular hexagonal carbon structure. Graphite is the 3-d structure. 

■^polycyclic aromatic hydrocarbon 
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1.7. Investigation 

In the present work, we investigate a subset of extreme carbon stars, ttiose wtiicli ex- 
hibit the ll/xm absorption feature. Through radiative transfer modeling, we investigate the 
nature of these dust shells. We use theoretical models and meteoritic data to limit the pa- 
rameter space and thus reduce the degeneracy within the model results. In addition, we look 
for correlations between observed parameters, such as those that define the ll/zm feature 
(strength, position, etc) as well as mass-loss rates and expansion velocities associated with 
the dust shells. Finally we determine timescales associated with the dust shells. 



2. Observations and Data Processing 



We investigated 10 extreme carbon stars observed spectroscopically by the ISO SWS all 
of which show evidence for an ~ ll/im absorption feature (see Table[2]and Figs. [2| [3H and [5|) . 
These s ources were chosen by searching the ISO archive for spectra of extreme carbon stars 
listed in lVolk et al.l (119921 ) and selecting those with an apparent 11/im absorpti on feature. In 



additi on, we used the color-color classification of "very cold" carbon stars by lOmont et al. 



( 119931 ) to identify further potential sources. Unfortunately most of the potential sources 
found in the color-color space (e.g. IRAS 17583-2291) were not observed by ISO SWS, and 
the IRAS LRS spectra are too low resolution and/or too noisy to be used in the present 
study. Four of our sources (IRAS 02408+5458, IRAS 19548+3035, IRAS 21318+5631, and 
IRAS 23166 +1655) were previously studied using ground-based obser vations and were found 
to be consistent with a self-absorbed SiC feature (jSpeck et al.l 119971 ). Two of these sources 
needed an extra absorption component at ~10/im, which were attributed to interstellar 
absorption. As discussed in ^ 11.41 these have since be en the source of some controversy 



(Clement et al. 



2005 



Pitman. Speck &: Hofmeisterll2006l ). Following the modeling efforts of 



Volk et all (120001 ) . ISpeck et al.l (120051 ) identified IRAS 00210+6221, IRAS 06582 4-1507, and 
IRAS 17534— 3030 as exhibiting an SiC absorption feature. The modeling of IVolk et al. 



( I2OOOI ) did not include SiC and did not attempt to fit the ll/xm absorption. Thus, division 
of the observed spectra by their respective RT model spectra revealed the llfim absorption 
feature. In this paper we present the discovery of three more potential SiC absorption 
features in ISO SWS spectra (IRAS 01144+6658, IRAS 03313+6058, and IRAS 22303+5950). 
These were discovered by searching the ISO archive for any extreme or "very cold" carbon 
stars as determined by their location in the IRAS color-color space. Those sources that fall 
within region III without OH maser e mission were examined. The l ocations of the sta rs in 
our sample in IRAS color-color space (Ivan der Veen fc Habinglll988l : lOmont et al.lll993[ ) are 
plotted in Figdl Interestingly, all sources except IRAS 19548+3035 plot along a line parallel 
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to and between the blackbody emission and the B{T, A,T) x A ^ emission hnes. 

The raw ISO data were extracted from the ISO data archive, and we used the Off- 
Line Processing (OLP) pipeline, version 10.1. Individual spectral sub-bands were cleaned 
of glitches (caused by cosmic ray particles) and other bad data sections. Next, they were 
flat-fielded, sigma-clipped (using the default values a = 3) and rebinned to the final spectral 
resolution {R = AA/A), which ranged fr om 200 to 700, dep ending on the scanning speed of 



the SWS grating during the observation (ILeech et al.ll2003l ). The final spectra are presented 
in Figure [21 which also shows the best-fittin 



blackbody continuum for each spectrum and 
the resulting continuum-divided spectra. The blackbody temperatures of the continua at 
listed in Table [31 The continuum-divided spectra clearly show an absorption feature in the 
10-13/im region, the basic parameters of which (barycentric position, peak-to-continuum 
ratio, full width half maximum; FWHM) are listed in Table [31 The excellent match between 
the overall spectrum and a single temperature blackbody suggests that we are seeing an 
isothermal surface within the dust shell. This represents the depth at which the shell becomes 
optically thick. The lack of extra emission at longer wavelengths suggests that any outlying 
dust is low enough in density to have an insignificant contribution to the overall emission. 

In addition to the ISO SWS spectra. Fig [3l [H and [51 shows the IRAS LRS spectra and 
the IRAS 12, 25, 60 and 100/xm photometry measurements. The difference in the flux levels 
between the IRAS and ISO data is not unexpected, since these stars are variable. However, 
the shape of the spectrum does not change significantly between the two observations, sug- 
gesting that changes in the stellar luminosity do not significantly impact the structure and 
composition of the dust shells. 

In order to determine the cause of the 10-13/im feature and the factors that govern 
its strength, width and position we have tabulated the barycentric position, feature-to- 
continuum ratio and equivalent width of the feature (Table [3]). In addition, we have also 
tabulated where the barycenter of the SiC feature would be if the short wavelength side of 
the absorption is due to silicate (as has been postulated, see § 11.41) . This, along with the 
feature-to-continuum ratio measured at 9.7 and 11.3/im, can be found in Table [31 



2.1. The "30/im" feature 

Another prominent spectral feature exhibited by our sample of extreme carbon stars 
is the so-called "30yum" feature. This feature is relatively common amongst carbon-rich 



"'Best fits are achieved by eye and proceed by examination of the continuum-divided spectra. 
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AGB stars, PPNe and PNe and was first discovered in Kuiper Airborne Observatory (KAO) 



spectra of carbon stars and PNe (IForrest et al. 1 198 11 
magnesium sulfide (MgS: e.g. ICoebel fc Moseley IQSsI : 



It has been widely attributed to 



Nuth et al. 1985: Omont et al. 1995 



Begemann et al.lll994j : iHony et al.ll2002l ). Modeling this feature is beyond the scope of the 



present work, but will be investigated in follow-up modeling. Our models make no attempt 
to fit the "30/im" feature. 



3. Radiative Transfer Modeling 



Radiative t ransfer modeling has been particularly useful in investi gating extreme carbon 
stars (see § Ol IVolk et al.lEggj : lGroenewegenlll99 5' 'Vol k et al.lboooh. We us ed the 1-D ra- 



diative transfer program DUSTY (jivezic fc Elitzuriil995i : iNenkova et al.ll200Cll ). to determine 
the effect of dust shell parameters on the emerging spectra from carbon stars. In all cases , 
the central star was assu med to be at 3000 K (typical for an AGB star). ISpeck et al.l (120001 ) 
and iDePew et al.l (120061 ) showed that changing this tempe rature by ±1000 K did not signif- 
icantly change the radiative transfer model's spectra (c.f. IVolk et al.ll2000l ). For simplicity, 
dust grains are assumed to be spherical. While DUSTY can include other grain shapes, this 
expands parameter space to create more degeneracy between models and is beyond the scope 
of the present work. 



3.1. Radial dust density distribution 



We assume a radial dust density distribution of 1/r^ which would reflect a constant mass- 
loss rate. This choice is somewhat controversial and certainly needs justifying. Current dust 
formation models suggest that carbon star mass- loss is expecte d to be modulated on several 
timescale s, especially th at of the pulsation cycle (IWoitkdl2006l . and references therein). Fur- 
thermore IWoitkd (120061 ) has suggested that the dynamics in the dust-forming zones around 
carbon stars lead to inhomogeneous dust formation, producing fine scale structure in the 
density of the dust enve lope. While puls ation shocks are predicted to have a strong effect 



on local conditions (e.g. ICherchneffI 120061) , this is not reflected in temporal changes in the 



IR spectra of carbon stars (jCorman et al.ll2008l ). Pre vious models of extrerne carbon stars 
have included steeper a density drop-off (see Table (U IVolk et al.l Il992l . l2000l : iGroenewegen 
19951 : IGroenewegen et al.lll998l ). which is meant to represe nt the increasi i ig mas s-loss rate 
associated with the onset of the superwind phase. However, IVillaver et al. J2002al lb[) showed 
that the hydrodynamics in the circumstellar she ll wipe out density structu r e and leave a 
shell with a 1/r^ density distribution. In addition Rowan- Robinson &: Harris ( )l983l ) showed 
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that carbon star spectra can be well-fitted using such a density distribution. 

In our models we assume that modulations in density have been wiped out or are 
unimportant in determining the spectrum at these high optical depths, as we are clearly 
seeing an outer dust shell surface. We also assume we are only detecting the dust that has 
formed since the onset of the superwind. We assume that there was a sudden increase in mass 
loss at some time in the last ~10000 years and that this mass-loss rate is now approximately 
constant, with small scale fluctuations being unimportant for dust properties. Consequently, 
the 1/r^ density distribution suffices. As will be seen in §IU the impact of assuming a steeper 
the density distribution is to remove dependence on shell size, and thus remove the ability to 
out limits on timescales. If we were to adopt a 1/r^ density distribution, all the dust would 
effectively be contained close to the star and would reflect the total mass lost over only a 
relatively short period. 



3.2. Modeling grain-size distributions 



The issue of choosing a grain size distribution is interesting and certainly requires 
more discussion. Our initial modeling studies used an MRN distribution (i.e., n{a) pro- 
portional to a~^, where n is the number of the grains in the si ze interval (a, a + da) and 
q = 3.5;amin = 0.005/im; and a^ax = 0.25/xm: iMathis et. allll977l ). This was chosen because 
as will become evident below, we do not actually know the grain size distribution and MRN 
is as plausible as any other. However, the MRN distribution was developed for interstellar 
dust where th e balance of formation and destru ction is different from AGB circumstellar 
environments. iDominik. Sedlmayr. fc Gaill (119891 ) suggested that the grain-size distribution 
create d in the circumstellar environ ments of AGB stars has a steeper power law (i.e. a^). 



while iKim. Martin &: Hendryl (jl994j . KMH) modified the MRN distribution to include an 
exponential fall-off term. 

We should also consider the "observational" evidence for range and distribution of grain 
sizes in carbon star outflows. Meteoritic presolar SiC grains from carbon stars have a 
huge grain-size distribution, ranging for 1.5 nm up to 26/im, with the majority of grains 



70%) falling in the .3-0.7/im range (see iDaulton et al.l l2003l : iBernatowicz et al.l 12006 



Clayton fc Nittlen |2004| ). Half the ma ss of the presolar SiC found in the Murchison me- 
teorite is in grains larger than 0.6/im ( Virag et al. 1992i). Car bon presolar grains follow a 
similar grain size distribution to SiC fiBernatowicz et al.ll2006l ). However, the sample may 
be biased in favor of large grains, which may be more apt to survive the journey through 
the ISM and incorporation into a solar system body. 
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Conventional theories of grain growth cannot produce the largest grains. Since AGB 
stars typically have mass-loss rates < 10~^ Mq yr~^ there should be an upper limit to the 
grain sizes of ~ O.l/xm. However, this assumes an isotropic distribution of material. The 
existence of meteoritic titanium carbide (TiC) within presolar carbon grains necessitates 
density inhomogeneities in the gas outflows (c.f., inhomoge nities caused by turbulence in 
hydrodynamic models of carbon-rich dust formation regions; IWoitkell2006l ). In addition, the 
grains must be ~ l/xm to produce the ll/zm feature, implying that there is a large population 
of small grains (b ut not necessarily p reclud ing large grains). Constraints on grain size were 
also discussed by iMartin &: Rogers! (119871 ) who found an upper limit to the grain size of 
~0.1/im based on polarizatio n meas urements of the famous carbon star, IRC+10216. On 
the other hand, iGroenewegenI (119971 ) used polarization measurements to limit the grain size 
to 0.1 < a < 0.35/im and suggested that very small grains (<80 nm) may not exist a round 
carbon stars (however, smaller SiC grains are found in meteorites). Meanwhile. I Jural (119941 ) 
argued the case for grains larger than l/xm in the circumstellar shells of IRC-l-10216. 

Towards the end of the AGB the onset of the supe rwind may lead to mass-loss rates 
as high as a fe w xlO~'^ M^d yr~^ (e.g. Hony et al.ll2003l ) which could translate into larger 



grains (see e.g. iBernatowicz et al.ll2006l . and references therein). However, the relationship 
between the evolution of carbon star s and the co n seque nt evolution of grain sizes in their 
circumstellar shells was discussed by ISpeck et al.l (120051 ) . They argued that the increased 
mass-loss rates at the end of the AGB phase lead to smaller, rather than larger grains as 
suggested by meteoritic evidence. Even the highest observed mass-loss rates cannot account 
for formation of titanium carbide (TiC) grains unless t he distribution of material is not 
spherically symmetric and density enhancements exist (IBernatowicz et al.l 120061 ). Such a 



distribution of material makes the concept of grain size distributions even more complex. 

The range of grain sizes used i n previ o us ra diative transfer modeling attempts varies 

hugely. The IR C-l-10216 models of iGriffinI (Il990l) required grain sizes limited to 5 < a < 

50nm . Similarly iBagnulo. Doyle fc Griffin Jl995h required grains smaller than 50nm. IGroenewegenI 
( 119971 ) also modeled IRC-l-10216 and found that the spectrum was best fit using a single grain 
size of 0.16/im ( rather than a dist ribution of grain sizes). The single-grain size approach was 
also adopted by lVolk et al.l (120001 ) who assume a single grain size of 0.1/xm. This is supported 
by the success of the early carbon star models of iRowan- Robinson fc HarrisI (119831 ) who also 
used a single grain size (0.1/im). However, since extreme carbon stars are expected to be the 
direct precursors for post-AGB objects, it may be more appropriate to consider the models 
of post-AGB stars. Such mod eling efforts have found they need gra in s up to millimetre- or 
even centimetre- in size (e.g. lJura et al.l l2000l : iMeixner et al.ll2002l ). iMeixner et al.l (120021 ) 
used a KMH-like distribution with an effective maximura grain size of 200/im (hereafter re- 
ferred to as KMH200). Meanwhile, ISzczerba et al.l (119971 ) modeled a carbon-rich post-AGB 
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object with a power-law distribution of grains in the range 0.01-1/im. 



Given the ran ge of possible g r ain si zes and distributions, it is difficult to know how best 



to model the dust. iGroeneweged (Il995l ) argued that the specific grain-size is not important 
as long as the grains are small enough that the absorption and scattering properties are 
independent of grain size. However, given the arguments for a population of large grains, 
limiting grains to smaller than O.l/im is unrealistic. 

In order to determine the effect of the choice of grain sizes on the model spectra we 
have generated models using five additional grains size distributions: 1) MRN-like with 

. = 0.25/im; g = 5 (as suggested by 



Dominik. Sedlmavr. & Gail 


1989); 2) only 0.1 


Volk et al.ll2000: 


Rowan- Robinson & Harris 


1983 



Groenewegenlll995l : lMartin &: Roger ] |l987l ): 



3) the dominant meteoritic grain sizes, i.e. O.l-l/xm only; 4) the standard KMH distribution 
(i.e. n{a) oc a-'^e''/^°, amin = 0.005/im, ao = 0.2/im; q = 3.5); and finally 5) KMH200: KMH 
with ao = 200/im. The results of this study are shown in Fig. [6l which shows the effect of 
changing the grain size distribution while keeping all other parameters constant. 

If we examine the differences in the spectra generated by changing the grain-size distri- 
butions, this can be understood in terms of the absorption efficiency of the grains and breadth 
of the grain size distribution. Changing from our default grain size distribution (MRN) to 
a modified MRN-like distribution with a steeper power law drop off (g = 5) as suggested by 
Dominik. Sedlmayr. fc Gail! (119891 ) makes very little difference to the model spectral energy 
distribution (SED). Likewise, switching from MRN to KMH has little effect on the overall 
shape of the SED. That the MRN, q = 5 and the KMH models are so similar is because 
they are basically weighted towards the same small grains. While the weighting is different, 
the same-sized small grains dominate the SED. However, changin g the size distributio n to 
include larger dust grains as used in models of post-AGB stars (e.g. iMeixner et al.ll2002[ ) has 
an major effect. The SED shifts to peak at much shorter wavelengths. This is because of the 
reduction in the number of small grains in order to include larger grains. The proportion of 
larger grains is small, but the removal of the smaller grains makes it possible for the stellar 
photons to penetrate deeper into the dust shell and provides a large population of warmer 
grains, resulting in a warmer SED. In both the KMH200 and 0.1/im cases there is a lack of 
very small grains which absorb a lot of stellar photons and change the temperature distri- 
bution (i.e. after the first layer of dust the temperature is lower, but if there are no small 
grains the stellar photons penetrate further.) Since size distributions similar to KMH200 are 
typically associated with long-lived dust disks, this particular distribution is not considered 
further. 



For the meteoritic grain size distribution, the short wavelength side of the SED is similar 



- 15 - 



to the default (MRN) model, but now the SED is much narrower. This can be explained 
by the very narrow range of grain sizes. There are no very small grains that can be easily 
heated and thus the shorter wavelength emission seen in KMH200 does not occur, but the 
lack of small grains also allows deeper penetration of stellar photon leading to a narrower 
temperature distribution. 

This grain size study suggests that for most adopted grain-size distributions, the result- 
ing SEDs will be equivalent and we assume the MRN distribution as "generic" . However, the 
meteoritic grain-size distribution can narrow the overall SED. For this reason, our modeling 
efforts concentrate on the generic (MRN) and "meteoritic" grain size distributions, where 
the "meteoritic" is taken to be an extreme amongst the range of reasonable grain-size as- 
sumptions. The impact of choosing different grain size distributions and the implications of 
these differences for the potential errors in our models will be discussed in § HI However, es- 
sentially, most grain-size distributions will yield the same results except for "meteoritic" and 
KMH200. In both cases changes in optical depth, inner dust temperature and/or relative 
geometrical shell thickness can be manipulated to fit the spectrum. There is a degeneracy 
in model fits due to the relationship between these three parameters that will be discussed 
in the next section. As will be seen in the § HJ the grain size effects cannot be ignored. In 
one source the need for "meteoritic" grain-size distribution is clear. 



3.3. Parameter space investigated 



In addition to grain sizes and radial density distribution, the variables investigated 
with DUSTY are the inner dust shell temperature (Tinner), optical depth (specified at 10/im; 
TiOAtm), dust composition, and the geometrical thickness of the dust shell, ^ = i^out/R-in, where 
Rin and Rout are the inner and outer radii of the dus t shell, respectiv ely. The optical constants 
for the dust components came from lPegouriel (119881 ) : lHannerl (119981 ) and lDraine &: Led (119841 ) 
for SiC, amorphous carbon and graphite, respectively. In nearly all cases it was possible 
to generate more than one model to fit the spectra, consequently we also investigate this 
degeneracy in parameter space and look for realistic ways to restrict it. 



3.3.1. Degeneracies in radiative transfer modelling 



There is a clear degeneracy in the models because of the relationship between certain 
parameters, e.g. optical depth and geom etric shell thickness, or optical depth and the temper- 
ature at the inner edge of the dust shell. Ilvezic &: Elitzurl (119971 ) discussed these degeneracies 
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and the relationship between the different input parameters in radiative transfer models, but 
we need to understand these relationship if we are to understand what our models mean. 

The wavelength-dependent optical depth, r\ is defined by: 



dTx = p{r)Kxdr (1) 

where p{r) is the density of the absorbing/scattering particles (i.e. the dust grains) and kx 
is the wavelength dependent opacity of the assemblage of particles along the line of sight. 



/■Rout 

r\= I p{r)Kxdr 

where Rin is the inner dust shell radius and i?out is the outer dust shell radius. While 
Kx is dependent on the density distribution of the grains, maintaining the size, shape and 
composition (and crystal structure) of the grains means that kx will not change significantly. 
For simplicity we assume kx remains constant. In our models we assume the density of the 
dust shell drops off as 1/r^ from the central star. In addition, our models specify the relative 
geometrical thickness of the dust shell as ,^ = Ro^t/Rm- Thus we get: 



T\ = Kx 



(2) 



The value of Ri^ is set by the values we input for the star's effective temperature and the 
inner dust radius (or condensation temperature. Tinner)- If we know the temperature and 
luminosity of the star, we can use a T(r) oc r~5 temperature distribution to determine the 
relationship between Tinner and -Rin and substituting into Eq. [2] we get: 



(3) 

Here we assume that the grains are blackbodies because they are largely carbon. In- 
cluding the albedo would allow for a more accurate calculation, but this would depend on 
detailed dust parameters (like crystal structure), and the error incurred by our assumption 
is small (i.e. significantly less than an order of magnitude). 

Therefore, according the Eq. [3l for a dust shell with constant relative shell thickness 
the optical depth should increase with the square of the inner dust temperature. Alter- 
natively, if the inner dust temperature is fixed, then increasing the relative shell thickness 



167ro" ^2 ^ 1 



should decrease optical depth a little (as — 1)/,^). This latter effect becomes negligible for 
large geometric sizes. 



3.4. Determining the Dust Condensation Temperature 



3.4-1- Theoretical dust condensation models 

As seen in the previous section, the value of the dust condensation temperature is a 
source of degeneracy on radiative transfer modeling. In order to reduce this degeneracy we 
turn to dust condensation theory to determine a theoretical dust condensation temperature 
(i.e. Tinner) that is appropriate for our stars. 



As discussed in § II. 6^ many presolar grains can be attributed to carbon stars and 
are valuable resources for investigation of dust formation regions. For example, in preso- 
lar grains titanium carbide (TiC) is found in the center of carbo n (C) grains frorn AGB 



stars, but only one SiC grain has been found coated in carbon (IClayton fc Nittlen 12004 



Bernatowicz et al.ll2006l ). Consequently there have been many studies of the theoretical con- 
densation sequence in Galactic carbon stars in attempts to constrain the physical parameters 
of the dust condensation regions. These studies showed how the condensation sequence of C, 
TiC and SiC is dependen t on various parameters, mo st notably C/0 ratio an d gaspressurqj 



3_pr 

( iLodders fc Fegleyl Il995l : ISharp fc Wasserburd Il995l ) . ISharp fc Wasserburgl (Il995l ) argued 
that if carbon forms at a higher condensation temperature, closer to the star than SiC, there 
is a significant decrease in the amount of carbon available in the gas, and thus SiC and C 
do not form simultaneously, resulting in naked SiC grains. Therefore, for Galactic sources, 
the condensation sequence in the majority of carbon stars should be TiC — C — SiC, in 
order to produce the coated TiC grains and uncoated (naked) SiC grains seen in the mete- 
oritic presolar gra in s samples. Observational ev idence for naked SiC grains is discussed by 
Speck et al.l (120051 ). ISharp fc Wasserburgl (119951 ) argued that from kinetic and stellar model 
considerations, dust grains should form in the pressure range 2 xlO^^ < P < 4 x 10^^ bar^. 



Lodders &: Fegleyl (119951 ) also modeled the effect of C/0 and pressure on the conden- 



sation sequence in carbon stars, as well as the effect of s-process and nitrogen abundances. 
They also briefly discuss the effect of metallicity. The general trends in condensation tem- 



^Gas pressure is a measure of the mass- loss rate (M) convolved with the photospheric temperature (T*) 
and outflow velocity (uoxp)- 



6 



The expected range of jE^as pressures in the dust formation zone for 0-rich AGB stars in the LMC is 



10-'^ bars - P - 10"'' bars (|Diikstra et aliboosh 
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peratures are: (1) all condensation temperatures decrease as the gas pressure decreases; and 
(2) At C/0 > 1 the condensation temperature of graphite increases with C/0 (for a given 
pressure and otherwise constant composition). 

Figure [7] shows how the condensation temperature of carbon and SiC vary with the gas 
pressure in the dust condensation zone. For the range plotted carbon always forms first 
from a cooling gas. Very high pressures are required for SiC to form before carbon. For 
solar metallicity and C/0 = 1.05, SiC forms before carbon for P > 3.4 x 10~^ bars. As 
C/0 increases, the minimum pressure required to form SiC first increases. Above C/0 ~1.5 
carbon always forms before SiC. The exact C/0 ratio at which the carbon forms before SiC 
depends on pressure. Therefore, in order to account for observations of SiC features in the 
Galaxy, and the presolar grain record, we can restrict the P-C/0 space such that, for low 
C/0 the gas pressure must remain low, but for higher C/0 the pressure can be higher. This 
can be used to constrain the dust forming environment around Galactic carbon stars. 

Studies of carbon star spectra in the Galaxy and the Magellani c Clouds have l ead to 
different interpretations with respect to the condensation sequence. iLagadec et al.l (120071 ) 
argued for a sequence in which SiC forms before C in the Galaxy, whereas SiC and C form 
together in the Large Magellanic Cloud (LMC) and the sequence is reversed (C, then SiC) 
in the Small Magellanic Cloud (SMC). However, they also suggest that the change in the 
strength of the features is due to the lower number of Si atoms available for SiC formation. 
The proposed Gala ctic condensatioii seque nce is at odds with both the models and the 
meteoritic evidence. iLeisenring et al.l (120081 ) support the condensation sequence in which C 
forms before SiC for the Galactic carbon stars, while finding that the Magell anic Clouds 



tended to form SiC first, followed by simu ltaneous condensation of SiC and C. ISpeck et al. 



( I2OO6I ) used the iLodders fc Fegleyl (119951 ) model to explain an unusual LMC carbon star 
spectrum, which suggests that the condensation sequence is sensitive to both metallicity and 
mass-loss rate. 



3.4-2. P-T space in the condensation zone around extreme carbon stars 

In order to constrain the input parameters to our model we need to be able to determine 
the pressure-temperature space in the dust condensation zone. For a mass-losing star with 
a mass-loss rate M and an expansion velocity of fexp, the density p of the circumstellar shell 
at a radius r is given by: 

M 



Anr'^Ve^p 
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If we know the temperature and luminosity of the star and the composition of the out- 
flowing material we can combine this information with the Ideal Gas Law and a T(r) oc 1/ 
temperature distribution to determine the gas pressure at the condensation radius, where 
the condensation radius is the distance from the star where the gas has the condensation 
temperature (Tcond)- 

For simplicity, the solid and g as phases are assume to be a t the same temperature. While 
this is clearly a simplification (e.g. IChigai fc Yamamotdl2003l ). the temperature difference is 
small compared to the difference needed to significantly affect dust formation. We assume 
that most of the outflowing material is atomic hydrogen. In fact it will probably be a mixture 
of atomic and molecular hydrogen (H2) since H2 forms around 2000 K and the temperature 
in the outflow is decreasing from the stellar surface temperature of ~3000 K to the dust 
condensation temperature in the 1000-1800 K range. An entirely molecular hydrogen gas 
would halve the gas pressure compared to the the atomic gas. However, we also assume 
published CO outflow velocities, which reflect the speed of the outflowing material after 
radiation pressure acceleration. Adopting the pre-dust-formation outflow speed (~5km/s), 
would increase the pressure. Thus using published values for our sample stars' luminosities, 
mass-loss rates and expansion velocities, we can estimate where their dust condensation 
zones fall in P-T space. 



3.4-3. Comparison of P-T space for dust condensation models and sample stars 

Figure [7] shows how the P-T space for the dust condensation zone for our target stars 
compares to condensation models in P-T space. It is clear the pressure is never high enough 
for SiC to form before carbon. Since our sample stars are expected to be the carbon stars 
with the highest mass-loss rates, it implies that this NEVER happens in Galactic carbon 
stars. This agrees with the meteoritic evidence which suggests that SiC does not get coated 
in carbon and supported the idea that variations in the strength of the SiC feature are related 
to self-absorption. 

The dust condensation temperature is dependent on the pressure in the gas from which 
the dust forms. In addition, the P-T space occupied by the target sources suggests that 
carbon grains will form at temperatures ^1600K. If the C/0 is very high, then graphite 
could form as high as ~1800K. However, Galactic carbon stars for which the C/0 ratios 
have been measured show that it is in the range 1 to 1.8, with a mean C/0 ratio of ~1.15 
Jlambert et al.lll986l : loiofsson et aDll993al jbl) . The precise C/0 for our sample is not known. 



but even low C/0 stars yield Tcond ^1550K (see Fig. [7]). Previous models of carbon stars 
have assumed much lower inner dust temperatures. The pressures and temperatures in the 
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gas around these stars meet the criteria for forming carbon dust at Tcond ^ 1600 K, (precise 
temperature depends on the C/0 ratio). Therefore it should form at these high temperatures. 

Once dust starts to form, the radial temperature profile will change due to absorption 
of starlight by dust grains. Therefore, we use the comparison above only to determine the 
inner dust radius. Using radial temperature profiles from our models we show that the 
temperature drops significantly more rapidly than 1/v^ (see Fig. [7]), which suggests that 
the next condensate (SiC) forms fairly close to the inner dust radius and thus mitigates the 
problem with DUSTY that the grains are assumed to be co-spatial. Indeed the models show 
that the temperature in the dust shell drops to the SiC formation temperature at ~1.3Rin 
for all cases, which is small compared to the shell thickness, even for the thinnest shells. 



3.5. Constraining Dust Composition 



Since our stars are carbon-rich, we limit the models to only including carbonaceous 
species such as graphite, amorphous carbon and silicon carbide. The choice of carbon grains 
is equivocal. We do not know whether the grains composed mostly of carbon are glassy, 
poly-nano crystalline or well-ordered. The meteoritic presolar "graphite" grains suggest that 
circumstellar dust can contain either well-ordered graphite or poly-nanocrystalline-graphene 
grains. It should be noted, however, that even the disordered graphene sheets are considered 
to be more graphitic than amorphous. 

Amorphous carbon essentially consists of a mixture of sp^ (graphite-like) and sp^ (diamond- 
like) carbon bonds. On heating, sp^ bonds tend to convert to sp^ bonds, thus graphitiz- 
ing the amorphous carbon, but it will remai n dense like the amorphous phase (2.8g/cm^ 
Comelh et allllQSsI : ISaadal I2OO0I : iKeliresI Eoosf ) . At ~1300K the fraction of graphitic bonds 



is ~90%. Comparison of the gas pressure in the circumstellar outflows to graphite formation 
temperature in § 13.4.31 shows that graphite should be able to form at temperatures signifi- 
cantly above 1300 K. Even if solid state carbon forms as a chaotic solid, at these temperatures 
it will quickly anneal to a graphitic form. Consequently, we argue that graphite is the better 
choice of carbonaceous material for modeling extreme carbon star dust shells. 



The hmitations of the use of the iPegouria (119881 ) data are discussed in iPitman et al. 



( I2OO8I ). Clearly, this data cannot produce the broad absorption feature seen in the observa- 
tions. However, we can use the relative changes in composition from star to star as a guide 
to understanding why these stars have the features we see. 
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4. Radiative Transfer Modeling Results 

The results of the radiative transfer modeling can be seen in Fig. [HI [9] and [TOl The 
parameters used in each case can be seen in Table [51 For each source, we have produced 
models using both MRN (generic) and "meteoritic" (extreme) grain size distributions. 

4.1. The effect of grain-size distribution 

Table [5] shows the model input parameters for both grainsize distributions. In all cases 
switching from MRN to "meteoritic" grainsizes leads to the need for increased SiC component 
(compared to graphite), typically requiring a three- or fourfold increase in the SiC fraction. 

There are other general trends including the need for decreased optical depths (rio^m) 
and increased geometric shell thickness. However, these trends do not hold for all objects. 
In the cases of IRAS 00210+6221 and IRAS 06582+1507, the only difference in parameters 
between the different grain size models is the fraction of SiC. IRAS 01144+6658 is the only 
source for which optical depth was increased and shell thickness was decreased in the "me- 
teoritic" model. 



4.2. Dust Shell Thicknesses 

As discussed in § 13.3.11 there is a degeneracy between relative shell thickness, inner dust 
temperature and optical depth. The inner dust temperature variability has been limited to 
1550-1800 K by theoretical considerations (see § 13.4.11) . It is possible to tweak parameters 
such as shell thickness, inner dust temperature and optical depth and get almost identical 
models. Fig. [TT] shows two almost identical model spectra for IRAS 23166+1655 with dif- 
ferent inner dust temperatures, optical depths and shell thicknesses. However, reducing the 
inner dust temperature merely requires reduction of the optical depth. While increasingly 
geometrically large dust shells can be accommodated by decreasing optical depth (see Eq. [3]), 
at some point, this also breaks down, as it leads to a significant population of colder grains 
which emit too much at long wavelengths. In this way we can use the models to place to 
an upper limit to the shell thickness. Table [H lists the published expansion velocities for 
our extreme carbon stars. From the model parameters, the physical value for i?out can be 
calculated, which gives the physical size of the dust shell. Since the models give an upper 
limit to the shell thickness, this upper limit to i?out together with the expansion velocity 
(■^exp) was used to calculate the time since the outermost edge of the shell was ejected from 
the star. The resulting ages of the dust shells are listed in Table [6] 
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As can be seen in Fig. |6l it is possible to accommodate thick shells, if we assume a much 
lower Tinner (1000 K), and assume that the dust grain are composed of amorphous rather than 
graphitic carbon, but these parameter values are precluded by the theoretical constraints on 
Tinner and compositiou/crystal structure discussed above. Clearly, since our constraints are 
theoretical, they may change as hypotheses are refined. 

The timescales for increased mass- loss are model dependent. As can be seen in Ta- 
ble O using a "meteoritic" grain size distribution generally requires a geometrically thicker 
shell, as well as lower values for rio^m and higher percentages of SiC. However, timescales 
associated with our dust shells are always very short (less than a few thousand years) re- 
gardless of grain-size distribution. Moreover, our timescales are consistent with those de- 
rived for the "superwind" seen in post-AGB objects (e.g. 



Skinner et al. 1997: Meixner et al. 



20041). These timescales a re too short to be associated with the theoretical superwind (e.g. 
Vassiliadis fc WoodI Il993l ). which is expected to last up to ~10% of the duration of the 
thermally-pulsing AGB phase. However, if we compare the number of extreme carbon stars 
to the total number of carbon stars in the Galaxy, we find that extreme carbon stars con- 
stitute only ~0.1% of the total C-star population. The thermally pulsing AGB phase is 
expected to last ~ 10^ years, and the time a star spends as a C-stars is even shorter (e.g. 
Lagadec et al.l 120071 ). Consequently, we might expect the dust-obscured phase to only last 
~ 10^ years. This is consistent with t he model timescale s , alth ough the MRN-models still 
appear to have very short timescales. iLagadec fc Zijlstral ( 120081 ) suggested that the trigger 
for the superwind is a combination of luminosity and carbon abundance. Although the du- 
ration of the C-star phase maybe be ~ 3 x 10^ years, for much of this time, a C-star will be 
below the critical carbon abundance required to drive the superwind and obscure the star. 

How does the extreme carbon star phase fit into the broader C-star evolutio n? Many, if 
not most, carbon-rich post-AGB s tars showed r narked axisymrnetric m orphologies (IMeixner et al, 
I999I : IWaelkens k WatersI Il999l : ISahail booi ISoker k Sub"a3 boosh . The extreme carbo n 



stars are expected to be the direct precursors of these objects (e.g. ISkinner et al.l Il998l ) 



but presently show little evidence for axisymmetr}0. The cause and timing of this axisym- 
metric structure is not known, but is believed to occur at the very end of the AGB phase. 
Dijkstra k Speckl (120061 ) showed that significant axisymmetry is not expected to develop un- 
til the last few tens or hundreds of years of the superwind phase. The onset of axisymmetry 
also leads to an optically thicker toroid of dust, as the circumstellar shell becomes equatori- 
ally enhanced. We suggest that the very short timescales associated with our model results 
may indicate that the extreme AGB stars are in the process of developing axisymmetry, but 



^one would expect to see more near-IR emission in the spectrum of strongly axisymmetric objects due to 
scattering of starlight escaping in the bipolar axis direction. 
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that this has not developed to the point of allowing large amounts of NIR scattered light 
into the spectrum. 



4.3. Dust Shell Density Distribution 

In §[3] we argued the case for maintaining a dust density distribution that follows a 
law. However, it has been argued that the increasing stellar luminosity and the onset of the 
superwind phase should give rise to a steeper density drop-off. 

Following the same arguments as shown in § 13.3. giving rise to Eq. [3] we can derive 
an equation in which the exponent of the density power law is a variable. This gives us the 
relationship between the optical depth (ta), inner dust temperature (Tinner), relative shell 
thicks (^) and the exponent of the density power law (x): 

[sj^TL.^ (4) 



(x - i)e- 



where kx is the wavelength dependent opacity of the assemblage of particles along the line 
of sight; is the luminosity of the star; and a is the Stefan-Boltzmann constant. 

As an example of the effect this, we assume the exponent, x = 3 (which is the highest 
value in the previous models). In this case, the relative shell thickness dependence dwindles, 
and the optical depth is essential strongly dependent on the inner dust temperature only. 
Since the previous models have severely underestimated the in ner dust temperatu re, they 
have also underestimated the optical depth (which may explain lGroenewegenlll995l ). 



The lack of dependence on shell thickness at x > 2 leads to a situation where most of the 
dust is essentially confined to a region close to the star and the outer dust becomes negligible, 
and so this confinement to the inner region is effectively the same as our assumption that 
there was a sudden increase in mass loss in the recent past. However, we attempted to 
model our sources with a radial density distribution and found that we cannot match 
the shape of the SED without restricting the geometrical shell thickness. Using the 
radial density distribution with the MRN grain-size distribution and -Rout/-Rin ~20 produces 
a SED that has too much emission longwards of ~ 20/im. Therefore, even if the radial 
density distribution is indicative of a steep increase in mass-loss rate the shells still need to 
be geometrically thin, and the shell thickness and percentage of SiC needed are essentially 
the same as for the models (see Fig. [11] for an example.) If the "meteoritic" grain- 
size distribution is used in conjunction with a 1/r^ drop-off, the observed SED cannot be 
matched. The model SED becomes too narrow, and increasing the outer dust-shell radius 
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does not help. 



4.4. IRAS 17534-3030 

IRAS 17534—3030 is exceptional amongst the present sample of extreme carbon stars, 
and certainly requires more discussion. As can be seen in the flux-calibrated and continuum- 
divided spectra in Fig. [2l IRAS 17534—3030 does not exhibit the usual molecular absorption 
bands around 13.7/im (due the C2H2) and shortwards of ~ 8yum (due to both C2H2 and 
HCN), which are present in the other sample sources. Its 11-13/im feature is intermediate 
between the narrow feature exemplified by IRAS 23166+1655 and the broad feature exem- 
plified by IRAS 19548+3035, indicating that at least some of whatever substance causes the 
broadening is present around IRAS 17534—3030. One of the sug gested carriers for th e short 



wavelength broadening of the 10-13/im absorption feature is C3 (jZijlstra et al.ll2006l ). How- 
ever, C3 is expected to be photospheric, rather than circumstellar, which probably precludes 
its detection in optica lly obscured stars. Moreover, the theoretical spectrum of C3 from 



J0rgensen et al.l (l2000l ) shows a strong absorption close to the ~ 5yum CO line, which is 
stronger than the ~ ll/im feature. The spectrum of IRAS 17534—3030 does not show this 
5/im absorption band. This, together with the lack of other molecular absorption feature 
provides evidence that the broadening of the 10-13/im feature is not molecular in origin. 

In addition to the lack of molecular absorption in its spectrum, IRAS 17534—3030 is 
unique is another way: it cannot be modeled with the MRN grain-size distribution. Mod- 
eling of this source requires the "meteoritic" grain-size distribution because of its narrow 
SED. Whereas our other sources can be fitted with either MRN or "meteoritic" grains, 
IRAS 17534-3030 cannot. 



4.5. Impact of the dust condensation temperature 

The increase in the inner dust temperature from ^ 1000 K in previous models up to 
^ 1600 K decreases the inner dust radius significantly. Consequently the flux of energy from 
the star hitting the inner dust radius is increased, leading to a greater acceleration and conse- 
quently a more effectively dust-driven wind. If graphite is formed by annealing of amorphous 
carbon, the grain density should remain high (with p=2.8g/cm^) and the increased acceler- 
ation is entirely due to increases flux of stellar photons. However, if graphite forms directly, 
rather than by annealing of amorphous carbon, the radiation pressure effect is further en- 
hanced, because graphite grains generally have a lower density than amorphous carbon. 
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Assuming 0.1/im-sized grains, the acceleration felt by graphite grains (with p=2.2g/cm^; 
c.f. meteoritic presolar grains; see § 11.6.21) at Tinner = 1800 K is 13 times greater than the 
acceleration felt by amorphous carbon grains (with p=2.8g/cm^) at Tanner = 1000 K. 



4.6. Correlations between observed and model parameters 

In the course of this investigation we have compiled a large number of parameters 
for these stars. For instance, Table [3] lists the peak position, peak strength, FWHM and 
equivalent width of the ~ 11/im absorption feature. Table H] lists the published mass-loss 
rates, luminosities and expansion velocities for our sample stars. 

In order to understand the physical conditions that give rise to the observed llpm ab- 
sorption feature, we have sought correlations between various parameters associated with 
the sample stars (as found in Table [21 IH E]). We have looked for correlations between each 
of the following: mass-loss rate (from CO observations); expansion velocity (from CO obser- 
vations); stellar luminosity; strength of the observed absorption; FWHM of the absorption; 
barycentric position of the absorption; equivalent width of the absorption; and model pa- 
rameters. 



Loup et al. J1993I ) showed that, for low mass loss rates there is a simple relationship 



between the [25]- [12] color and mass- loss rates. This breaks down at mass-loss rates ^ 
lO~^M0yr~^. It has been suggested that such high mass-loss rate stars at CO-emission 
deficient due to either saturation effects, low kinetic temperatures or possibly dramatic recent 
increases in mass-loss. If we extrapolate the trends from low mass-loss rates to determine 
mass-loss rates from the [25]-[12] color we find that these objects should have a mass-loss rate 
in excess of lO~^M0yr~^, consistent with the high modeled optical depths. Because of this 
relationship, we sought correlations between the various observed and modeled parameters 
and the [25]- [12] and [60]- [25] colors. This search yielded only one correlation: between the 
[60]- [25] colors and the dust mass- loss rate from previous models. It is possible t hat a better 



corre l ation may be found using the Manchester Method i.e. the [6]- [9] color (e.g. ISloan et al. 



20061 : IZijlstra et al.l 120061 ) . however, we suspect that the lack of correlation arises because of 



the intrinsic degeneracy in the modeling. 

In § 11.41 and §[2]we discussed the possibility that the broadening of the 10-13/xm feature 
might be due to silicate dust. With this in mind, we re-measured the position and strength 
of the llpm feature assuming that the short wavelength wing is due to silicate. This involved 
measuring the feature-to-continuum strength at 9.7 and 11.3/im. These data are tabulated 
in Table [3] and were also included in the investigation of correlations between parameters. 
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One pair of parameters that yielded a correlation were the best-fit blackbody tempera- 
ture and the model optical depth, which is demonstrated in Fig [121 This correlation occurs 
whether we assume MRN or "meteoritic" grainsizes. Since there is a relationship between 
the optical depth and shell thickness, this correlation seems to support models validity. 

In addition to this relationship between the modeled optical depth and the best-fit 
blackbody temperature, we found two other parameters that the optical depth correlates 
with: modeled percentage of SiC in the dust shell; and the calculated timescale of the 
obscuring dust. In both cases, these correlations only hold for the generic (MRN) grain- 
sze distribution models, (shown in Fig. [T3|) . The correlation between optical depth and 
percentage SiC is such that lower optical depths require more SiC. This in turn suggests 
that as mass-loss rates increase, the SiC component dwindles. This can be interpreted in 
two ways: (1) at these high mass-loss rates SiC gets coated by carbon; and (2) increased 
mass-loss rates are associated with higher C/0 ratios; carbon is enriched but not silicon, 
and thus more carbon grains can be made, but not more SiC grains. The first option seems 
unlikely in light of meteoritic evidence (§ 11.61) and theoretical condensation models (§ I3.4.ip . 
However, it is possible that there is a metallicity effect in play. The presolar grains were 
formed prior to the formation of the so l ar system, and thus the ir source stars may have had 



lower metallicities. ISpeck et al.l (120061 ): iLeisenring et al.l (120081 ) argued that coating of SiC 



grains is more likely in lower metallicity environments, and thus it is difficult to see how 
the highe r metallicity obiects we are now witnessing could have carbon-coated SiC grains. 



However, ILeisenring et al.l (120081 ) argued that SiC grains form the nucleation seed for MgS, 
which may explain this correlation. The second option may be logical if outflows are dust 
driven. If more carbon grains can form, the radiation pressure driving the outflow is more 
effective. 

The second optical depth correlation is with the calculated timescale (also shown in 
Fig. [island again, only for the MRN grain-size distribution). In this case, the timescale 
for the dustshell increases with optical depth. The oldest shells have the highest optical 
depth. This may simply be due to the pathlength dependence of optical depth (see Eqn. [T|), 
since r depends on the geometric size of the dust shell the oldest shells should have the 
largest optical depth. However, this should also lead to a correlation between model dust 
shell thickness and optical depth. No such correlation exists. Another interpretation is that 
older stars have higher mass-loss rates, and thus denser dust shells and higher optical depths. 
However, this assumes the stars all had the same initial mass. This apparent correlation is 
not strong, and does not holds for "meteoritic" grains. Consequently, we will not attempt 
to place too much importance on it. 

There is no correlation between the strength of the 9.7/xm absorption and that at 11.3/xm, 
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implying that the carriers of these features may not be related. Furthermore, there is no 
correlation between the 9.7/im strength and distance to the object, which leads us to suggest 
that the cause is not interstellar. Obviously there is a strong correlation between the strength 
of the 9.7/im absorption and the barycentric position of the feature, indicating that this is 
the cause of the broadening. 

The most interesting correlations are between the "isolated" 11.3/zm feature (i.e. the 
residual feature after taking silicate absorption into account) and the fitted blackbody tem- 
perature. There is a correlation both between the SiC barycentric position and the blackbody 
temperature and between the 11.3/im feature-to-continuum ratio (i.e. SiC feature strength) 
and the blackbody temperature (both shown in Fig. [T2|) . 



The position of the SiC absorption feature moves to shorter wavelengths for higher 
blackbody temperatures while also becoming weaker. This can be understood as being the 
result of more absorption occuring when the surface we see is at a lower temperature, which 
is associated with longer wavelength absorption. 



We can interpret this result in terms of the self-absorption scheme described by lSpeck et al 



(I2OO5I ). The shifts in wavelength were attributed to a changes in grain size. However, the 



weaker absorption and warmer blackbody temperature s associated with t he shortest wave- 



length peaks is the opposite trend to that described by ISpeck et al.l (120051 ). Following their 
scheme, this would imply that the weakest absorption is associated with the largest SiC grain. 
The problem with interpreting this observation is that we do not know the C/0 ratios for 
these stars. While C/0 does not affect the SiC condensation temperature, it will affect the 
graphite condensation temperature. As seen in Fig [111 ^6 can reduce the inner dust tem- 
perature and compensate by decreasing both the optical depth and the shell thickness. This 
would also change the depth into the shell at which the SiC forms. It is also possible that 
the shift in position is due to incomplete subtraction of the silicate contribution. Further 
studies are ongoing, but are beyond the scope of the present work. 

The assumption that there is a silicate contribution to the spectrum raises several ques- 
tions. Identifying the source of the 9.7/im contribution to the absorption feature remains 
beyond the scope of the present paper, but, needless to say, if it is silicate, the specter of dual 
chemistry wo uld require seriousl y reth inking our current models of dust formation around 



carbon stars. iLodders fc Fegleyl (119991 ) suggested that silicates could form around carbon 
stars, and indeed silicate carbon stars exist. It is possible that the feature may be attributable 
to some form of hydrogenated amorphous carbon, which has been p ostulated as the source 



of an ~ 9/im emission feature in optical thin carbon star spectra (see [Thompson et al.ll2006 
and reference therein). 
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No additional correlations were found. The lack of correlation be t ween most of the 
observable or modeled parameters echoes the results of [Thompson et al.l (120061 ) . who found 
that there are no trends in the parameters associated with ll/xm emission feature in vis- 
ibly observable carbon stars. This suggests that even amongst the extreme carbon stars, 
variations in C/0, s-process and nitrogen enhancements make discernment of the physical 
properties associated with the 11/xm feature very difficult. 



5. Conclusions 

We have presented three previously unrecognized SiC absorption features in the spectra 
of extreme carbon stars. Together with the seven known SiC absorption stars, this brings 
the total of known extreme carbon stars with SiC absorption features to ten. 

Previous radiative transfer models of extreme carbon stars utilized relatively low con- 
densation temperatures. Here, theoretical condensation models have been used to justify 
much higher condensation temperatures. In addition, our models use graphite instead of 
amorphous carbon, because of the preferential formation of graphite at higher temperature 
and the meteoritic evidence. Both the higher condensation temperature (through a decrease 
in the inner radius of the dust shell) and to a smaller extent the use of graphite will greatly 
increase the acceleration felt by the dust grains in the shell relative to parameters used in 
previous research. 

We have shown that grain-size issues cannot be ignored in the production of models 
that accurately fit the observed spectra of extreme carbon stars. The size distribution that 
is needed is not clearly defined because of the inherent degeneracy in radiative transfer 
modeling. Meteoritic grain-size distributions are as valid as other size distributions with the 
advantage of being model independent. However, they may underestimate the contribution 
from small grains. 

With the exception of IRAS 17534—3030, all sample stars could be modeled with either 
the generic (MRN) or "meteoritic grain-size distribution IRAS 17534— 3030's narrow SED 
required the use of the "meteoric" distribution. Furthermore, there is no evidence for any 
molecular absorption in its spectrum. Because the ll/xm feature is still present in the absence 
of the other molecular features, it supports the attribution of this feature to a solid state 
carrier. 

The various parameters compiled in the course of this research (both through radiative 
transfer modeling and from observations) have been compared in order to identify any cor- 
relations, with the result that the cause of differences is the spectra cannot be attributed to 
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mass-loss rate or gas pressure in the dust conde nsation zone. In f act th e paucity of corre- 



lations between parameters echoes the results of [Thompson et al.l (120061 ) and suggests that 
even amongst the extreme carbon stars, variations in C/0, s-process and nitrogen enhance- 
ments make ll/im a poor probe of the details of dust shell parameters. 

The timescales associated with the heavy mass-loss experienced by these extreme carbon 
stars are very short (tens to hundred of years) and are not consistent with timescales for 
the superwind. This indicates that the heavy mass-loss phase of carbon stars is not a direct 
result of thermal pulse (although thermal pulses may be the root cause). 

This work is supported by NSF AST-0607341. We are very grateful to the referee, 
Albert Zijlstra, for his comments which have significantly improved this paper. Kevin Volk 
is also thanked for helpful advice on this paper. 
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Fig. 1. — Regions in the IRAS color-color diagram populated by extreme carbon stars, see 
text for details. The black body curve is indicated by the solid grey line. The modified 
blackbody [-6(7", A) * A~^] is indicated by the dashed grey fine 
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Fig. 2. — ISO SWS spectra of ten extreme C-stars. Left Panel: Flux-calibrated spectra and 
fitted blackbody continua. Solid line = ISO spectrum; dashed line = best-fitting blackbody 
continuum; The y-axis is the flux (F^) in Wm~^ fim~^] the x-axis is wavelength in fim. Right 
Panel: Continuum-divided spectra. Table [3] lists the blackbody temperatures used in each 
case to produce the continuum and the continuum-divided spectra. 
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Fig. 3. — ISO SWS spectra of sample extreme C-stars together with IRAS 12, 25, 60 and 
100/im photometry points and IRAS LRS spectra. SoUd hne = IRAS spectrum; points = 
IRAS photometry points; dashed hne = ISO SWS spectrum. The |/-axis is the flux (AF^) 
in Wm~^; the x-axis is wavelength in /im. 
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Fig. 4. — ISO SWS spectra of sample extreme C-stars together with IRAS 12, 25, 60 and 
100/im photometry points and IRAS LRS spectra (part 2). Sohd hne = IRAS spectrum; 
points = IRAS photometry points; dashed line = ISO SWS spectrum. The y-axis is the 
flux (AF;s^) in Wm~^; the x-axis is wavelength in /xm. 
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Fig. 5. — ISO SWS spectra of sample extreme C-stars together with IRAS 12, 25, 60 and 
100/xm photometry points and IRAS LRS spectra (part 3). Sohd Une = IRAS spectrum; 
points = IRAS photometry points; dashed hne = ISO SWS spectrum. The y-axis is the 
flux (AFa) in Wm-2; the a;-axis is wavelength in jj,m. 
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Fig. 6. — The effect of grain-size distributions on tlie model spectral energy distribution. In 
all cases the model parameters are identical except for the grain-size distribution, see § 13.21 
for detailed description of the grain-size distributions. The grey line shows the ISO-SWS 
spectrum of IRAS 03313+6058 for comparison. 
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Fig. 7. — Pressure- Temperature space for the dust condensation region around extreme 
carbon stars (assumes solar metallic ity). Grey lines indicate the condensation temperature 
for a given pressure as calculated by iLodders fc Fegleyl (119951 ) . Black lines indicate the P-T 
paths for the otuflowing gas from our sample stars; (top) as calculated from the published CO 
mass-loss rates and expansion velocities (see § 13.4.21 for details); and {bottom) as calculated 
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Fig. 8. — Best fit models (part 1). solid line = ISO-SWS spectrum; dashed line = best fit 
model using MRN grain-size distribution; dotted line = best fit model using "meteoritic" 
grain-size distribution; X-axis is wavelength (yum); y-axis is fiux (AFx) in Wm~^. In all 
cases, Tj,=3000K and Tinner=1800 K. Model parameters are listed in Table [5] 
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Fig. 9. — Best fit models (part 2) solid line = ISO-SWS spectrum; dashed line = best fit 
model X-axis is wavelength (/im); y-axis is fiux (AFx) in Wm~^. In all cases, T^=3000K and 
^inner=1800 K. Modcl parameters are listed in Table [5] 
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Fig. 10. — Best fit models (part 3) solid line = ISO-SWS spectrum; dashed line = best fit 
model X-axis is wavelength (yum); y-axis is fiux (AFx) in Wm~^. In all cases, T^=3000K and 
^inner=1800 K. other parameters are indicated in the legend. Parameters of these models are 
also listed in Table [5] 
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Fig. 11.— Degeneracy in the best fit models for IRAS 23166+1655. solid line = ISO-SWS 
spectrum; dashed fine = best fit model with Tinner^lSOO K dotted line = best fit model with 
21nner=1600K X-axis is wavelength (/im); y-axis is fiux {XFx) in Wm~^. In both cases, 
T^=3000 K other parameters are indicated in the legend. 
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Fig. 12. — Correlations between the best- fit blackbody temperature and (top) modeled op- 
tical depth; (middle) 11.3//m feature-to-continuum ratio; and {bottom) SiC feature position. 
Solid lines represent the linear regression fit through the points. In the top figure, the dotted 
line represents the linear regression for model fits with MRN grainsize distributions. Grey 
points are for meteoritic grainsize models, and grey dashed line is the linear regression fit 
through the points. 
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Fig. 13. — Correlations between the modeled optical depth and (top) the modeled %age SiC 
( bottom) the calculated superwind timescale from Table [6] Solid lines represent the linear 
regression fit through the points. 
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Table 1: Previous Models of Extreme Carbon stars with ll/im absorption features 



Star 


composition^ 


Grain-size 


TiLSfim drop-off 


R*„ Tj„ M* 


REF 




(SiC%) 


(^m) 






(xlO-^pc) (K) Meyr-i 




00210+6221 





0.1 


4.5 


3.00 


0.362/D 1000 12x10" 


2 


01144+6658 


8 


0.1 


4.84 




1000 9.50 xlO-'^ 


4 


06582-1507 







5.15 




1000 1.08x10"^ 


1 







0.1 


2.1 


2.25 


1.23 1000 2.59x10"^ 


2 


17534-3030 





0.1 


4.4 


2.50 


1.39 1000 8.20x10"^ 


2 


19548+3035 





0.1 


2.5 


2.50 


1.17 1000 3.89x10-^ 


2 


21318+5631 







1.36 




700 1.10x10"^ 


3 


23166+1655 





0.1 


1.12 




650 5.50x10-^ 


4 




30 




7.85 




1000 3.30x10-5 


1 









1.19 




650 1.00x10-^ 


3 



References 1: IVolk et al l (119921 ): 2: IVolk et al.l (|2000l ): 3: iGroeneweeenl (|l995l ): 4: 
Groenewegen et al.l (jl998l ) 

* originally quoted as a fu nction of distance. Value quoted here assume distances from 
Groenewegen et al.l (j2002l ) , listed in Table H] 



Neither Tvolk et al.1 
models 



nor 



t composition assumes reniainde r dust is carbon. In all but lVolk et al.l (119921 ) the carbon 
is amor phous: IVolk et alJ (119921) uses graphi t ic car bon. 



Groenewegen! (|l995l ) specify the grains sizes used in their 



Table 2: Target List 
IRAS Other 
Number Names 



R.A. Z Decl. TDT Date of 

(J2000) (J2000) number Observation 

40401901 1996 Dec 24 
68800128 1997 Oct 03 
80002504 1998 Jan 24 
62301907 1997 Jul 31 
71002102 1997 Oct 26 
12102004 1996 Mar 17 
56100849 1997 May 30 
11101103 1997 Mar 7 
77900836 1998 Jan 02 
37900867 1996 Nov 29 



00210+6221 CGCS 6006 00 23 51.2 +62 38 16.4 
01144+6658 V829 Gas, AFGL 190, GGCS 6017 01 17 51.6 +67 13 55.4 

02408+5458 02 44 25.2 +55 11 15 

03313+6058 CGGS 6061 03 35 30.7 +61 08 47.2 

06582+1507 CGGS 6193 07 01 08.44 +15 03 39.8 

17534-3030 AFGL 5416, GGGS 6690 17 56 33.1 -30 30 47.1 

19548+3035 AFGL 2477, GGGS 6851 21 50 45.0 +53 15 28.0 

21318+5631 AFGL 5625S, GGGS 6888 21 33 22.98 +56 44 35.0 

22303+5950 GGGS 6906 22 32 12.8 +60 06 04.1 

23166+1655 LL Peg, AFGL 3068, GGGS 6913 23 19 12.39 +17 11 35.4 
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Table 3: Observed parameters of the target sources 
IRAS Feature Feature^ FWHM Equivalent dJJm? ll.S^m^ SiC^ 





(K) 


RarvppTit.pr 


ux v^xxc, uxx 




Width 

V V ±KX UXX 


strength 


btrengtn 


position 


00210+6221 


290K 


10.45 


0.73 


2.28 


0.60 


0.76 


0.77 


11.3 


01144+6658 


280K 


10.51 


0.82 


2.32 


0.35 


0.88 


0.87 


11.3 


02408+5458 


250K 


11.28 


0.73 


2.20 


0.30 




0.83 


11.3 


03313+6058 


350K 


10.12 


0.87 


1.72 


0.17 


0.90 


0.97 


10.8 


06582+1507 


320K 


10.25 


0.91 


0.70 


0.1 


0.97 


0.94 


11.0 


17534-3030 


280K 


10.59 


0.87 


0.80 


0.30 


0.92 


0.83 


11.3 


19548+3035 


295K 


10.24 


0.87 


2.39 


0.25 


0.88 


0.93 


11.0 


21318+5631 


300K 


10.20 


0.84 


2.90 


0.37 


0.88 


0.93 


11.0 


22303+5950 


350K 


10.25 


0.84 


1.70 


0.25 


0.86 


0.94 


11.0 


23166+1655 


290K 


11.42 


0.86 


2.00 


0.21 




0.88 


11.4 



Feature strength is the "peak"-to-continuum ratio and is measured and the barycentric position. 
^ The 9.7;um strength is the feature-to-continuum ratio measured at 9.7;um. 
^ The 11.3/im strength is the feature-to-continuum ratio measured at 11.3/im. 
^ The SiC position, the the approximate barycentric position that the SiC feature would have if 
the short wavelength side of the absorption were due to silicate 9.7;um absorption. 



Table 4: Compilation of CO mass-loss rates and expansion velocities, together with distances 
and luminosities of sample stars. 



Source Name 


^exp 


D 






Mdust 




L. 




km/s 


kpc 


Moyr-i 




Moyr-i 


Lo 


00210+6221* 


16.7 


3.97 


3.02x10" 


^5 






1x10^ 


01144+6658 


18 


2.78 


6.38x10- 


-5 


1.51x10" 


-7 


1.69x10^ 


02408+5458* 


11 


5.3 


1.60x10- 


-5 


3.5 xlO- 


7 


5.70x10^ 


03313+6058 


13.9 


5.24 


2.37x10- 


-5 


1.14x10" 


-7 


1.31x10^ 


06582+1507 


13.7 


4.7 


1.43x10- 


-5 


1.07x10" 


-7 


1.32x10^ 


17534-3030 


19 


2 


3.76x10" 


^5 


1.06x10- 


-7 


1.21x10^ 


19548+3035 


22.3 


3.38 


1.14x10" 


A 


2.15x10- 


-7 


1.32x10^ 


21318+5631 


19.6 


1.77 


7.69x10- 


-6 


1.18x10- 


-7 


1.24x10^ 


22303+5950 


18.3 


3.86 


3.19x10" 


A 


1.09x10- 


-7 


1.25x10^ 


23166+1655 


15.1 


1 


1.44x10" 


-5 


8.27x10" 


-8 


1.10x10^ 



are from iGroenewegen et al.l (120021) except those rnarked 



All data 

with *. IRAS 02408+5458 data come from IGroenewegen et al.l (Il999l ). 



IRAS 00210+6221 data are compiled fromlVolk et al l J200oh (M, whi c h also 
convolves the distance from Groenewegen et al. 2002 ): and lVolk et al. ( 1992 ) 



-'exp 7 
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Table 5: Model and Derived Parameters 



IRAS 


T. 


T 

-■- inner 


-^out/-^in 




Composition 


Number 










SiC% 


Graphite % 




MRN grain size distribution 




00210+6221 


3000 


1800 


10 


6 


10 


90 


01144+6658 


3000 


1800 


15 


6.5 


10 


90 


02408+5458 


3000 


1800 


20 


12 


3 


97 


03313+6058 


3000 


1800 


15 


4 


10 


90 


06582+1507 


3000 


1800 


10 


6 


5 


95 


17534-3030 


3000 


1800 










19548+3035 


3000 


1800 


15 


6.5 


5 


95 


21318+5631 


3000 


1800 


10 


8 


3 


97 


22303+5950 


3000 


1800 


15 


4 


10 


90 


23166+1655 


3000 


1800 


15 


7.5 


3 


97 




Meteoritic . 


grain size 


distribut 


ion 




00210+6221 


3000 


1800 


10 


6 


30 


70 


01144+6658 


3000 


1800 


10 


7.5 


30 


70 


02408+5458 


3000 


1800 


100 


8.5 


12 


88 


03313+6058 


3000 


1800 


500 


3.25 


40 


60 


06582+1507 


3000 


1800 


10 


6 


20 


80 


17534-3030 


3000 


1800 


20 


5 


35 


65 


19548+3035 


3000 


1800 


100* 


4.75 


25 


75 


21318+5631 


3000 


1800 


100 


5 


20 


80 


22303+5950 


3000 


1800 


500 


3.25 


40 


60 


23166+1655 


3000 


1800 


50 


5 


25 


75 



* Models with i?out/-Rin > 100 can be accommodated because the data beyond 26/um is poor 
and ignored, but based on the similarity of this source to IRAS 21318+5631, we assume this 
is a good upper limit. 
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Table 6: Size and age of the dust shells 

MRN size distribution Meteoritic size distribution 



IRAS 


Rin 


R-out 


Age 


Rin 


R-out 


Age 


Number 


(lO^^cni) 


(lO^^cni) 


(yrs) 


(lO^^cm) 


(lO^^cm) 


(yrs) 


00210 


2.68 


26.8 


50.9 


2.78 


2.78 


52.8 


01144 


2.74 


41.1 


72.4 


3.08 


30.8 


54.2 


02408 


4.11 


82.2 


236.8 


2.99 


299 


861.3 


03313 


2.17 


32.6 


74.2 


2.10 


1050 


2393.7 


06582 


2.73 


27.3 


63.1 


2.72 


27.2 


62.9 


17534 


2.51 


50.2 


83.7 


2.51 


50.2 


83.7 


19548 


2.78 


41.7 


59.3 


2.38 


238 


338.2 


21318 


3.18 


31.8 


51.4 


2.41 


241 


389.6 


22303 


2.24 


33.6 


58.2 


2.10 


1050 


1818.2 


23166 


2.88 


43.2 


90.7 


2.45 


123 


257.1 



Calculation of Rout is based on the model parameters listed in Table [51 Rjn comes 
from the model output. The calculation of the age of the dust shells is then done 
using these data and the observed expansion velocities hsted in Table H] 



